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Abstract 

These lecture notes give a very short introduction to coarsening phe- 
nomena and summarize some recent results in the field. They focus on 
three aspects: the super-universality hypothesis, the geometry of growing 
structures, and coarsening in the spiral kinetically constrained model. 



1 Introduction 

These notes are complementary to my lectures at FPSP XII. My aim is 
to discuss in this text some recent developments in the theory of coarsening 
phenomena [1] (although I only covered some of them in the talks). Coarsening 
is possibly the simplest example of macroscopic non-equilibrium relaxation and 
it has very far-reaching practical applications. This problem, although pretty 
old and rather well understood qualitatively, is still far from having a close 
and satisfactory quantitative description going beyond the pretty successful but 
somehow deceiving use of the dynamic scaling hypothesis. 

After presenting the models I shall focus on and the coarsening phenomenon, 
I shall discuss the following aspects: (1) super-universality; (2) the statistics 
and geometry of evolving ordered structures; (3) domain growth in a kinetically 
constrained spin model, the two dimensional spiral model. Some reasons why 
we studied these problems are the following. We wished to check the super- 
universality hypothesis in more detail than what it had been done so far. In our 
opinion, establishing the validity of this hypothesis might be useful to eventually 
develop a successful analytic method to compute scaling functions in coarsening 
problems. We analyzed the geometry of domain growth with similar ideas to 
the ones used in the study of critical equilibrium structures in the past. Know- 
ing in detail the time-dependent distribution of structures and their geometric 
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properties might give us hints into what to search for in more complex systems 
with out of equihbrium dynamics such as glasses. In this respect, kinetically 
constrained models are supposed to be toy models for glassy relaxation. We 
demonstrated that the mechanism of relaxation in the spiral model is of the 
coarsening type. This is interesting per se and for applications to glassy physics 
too. In short, apart from making our understanding of coarsening more com- 
plete, we expect that these results might be useful to better grasp the more 
complex dynamics of glassy systems. 

References to the original articles will be given in the main text. My results 
on this field have been obtained in collaboration with the colleagues and students 
that I warmly thank in the acknowledgements. 

2 Models 

In this section we introduce three kinds of problems with phase transitions: 
purely geometric models with percolation as the standard example; models with 
an energy (or cost) function that have both static and dynamic transitions, 
with Ising and Potts spin models as typical instances; and models with purely 
dynamic phase transition such as kinetically constrained spin systems focusing 
on the bi-dimensional spiral model. Related to the latter case, I shall also 
mention coupled map lattices as well as spin models with dynamic rules that 
do not satisfy detailed balance. 

2.1 Geometric and probabilistic models 

We recall the definition of two purely geometric and probabilistic models: 
site and directed percolation. These models change behaviour at a special value 
of their control parameter, pc, and this threshold has many points in common 
with usual thermal second-order phase transitions. Percolation theory deals 
mainly with the critical phenomenon and the study of cluster sizes and their 
geometric properties. 

2.1.1 Site percolation 

Percolation is the propagation of activities through connected space. The 
site percolation model [2] is defined as follows. For each site on the considered 
lattice one tosses a coin independently. With probability p the site is occupied 
and with probability 1 — p it remains empty. A configuration is constructed by 
sweeping the lattice once and filling the sites (or not) independently with this 
rule. 

Clusters are defined as groups of first-neighbour filled (or empty) sites. At 
a (lattice dependent) threshold Pc an infinite cluster (on an infinite lattice) 
appears. Strictly speaking, for p < Pc (1 < Pc) there is no infinite connected- 
component of filled (empty) sites with probability one while for p > pc {1 — 
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p > pc) there is an infinite connected-component of filled (empty) sites with 
probability one. The critical values of in some typical bi-dimensional lattices 
are: pc = 1/2 (triangular), pc « 0.59 (square), pc ~ 0.62 (honeycomb). 

Percolation can be obtained as a limit of the Potts model (sec its definition 
below) with q — > 1 [5] . 

2.1.2 Directed percolation 

Directed percolation [3] is percolation with a special direction along which 
the activity propagates. It mimics filtering of fluids through porous materials 
along a given direction. For example, in the percolation of coffee making, the 
source is the top surface of the grounded coffee that receives the water and 
the activity is to have water flow. On a square lattice the model is defined by 
assigning a directed bond on horizontal (say, pointing to the left) and vertical 
(say, pointing down) edges with probability p. This mimics different microscopic 
pore connectivity. This model displays a phase transition from a macroscopically 
permeable (percolating) to an impermeable (non-percolating) state. Directed 
percolation on a 2d square lattice has a transition at Pc — 0.705. 

More generally, the terms percolation and directed percolation stand for uni- 
versality classes of continuous phase transitions which arc characterized by the 
same type of collective behavior on large scales. A number of critical exponents, 
linked to the size of the percolation cluster, correlation functions, and so on can 
be defined similarly to what is done in thermal critical phenomena [2] . 

2.2 Hamiltonian systems 

Physical models, as well as some mathematical problems such as those re- 
quiring optimisation, have an energy or cost function associated to each con- 
figuration. In equilibrium, the configurations are sampled with a probability 
distribution function that depends on their energy and the external parameters 
such as temperature or a chemical potential. The experimental conditions, that 
is to say whether the system is isolated or in contact with thermal or particle 
reservoirs, dictate the statistical ensemble to be used (micro-canonical, canon- 
ical or gran-canonical). Out of equilibrium, say when the system is let evolve 
from an initial condition that is not one selected from the equilibrium measure, 
the system wanders in phase space in a manner that we shall discuss below. 

2.2.1 The Ising model 

Anisotropic magnets are modeled in a very simple way by using Ising spins. 
Si = ±, to describe the magnetic moments and by choosing an adequate inter- 
action among them. The Hamiltonian of the d dimensional Ising model is 
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The spin variables are placed on the vertexes of a lattice and the sum runs over 
nearest-neighbours. The exchange interactions Jij are all positive, favouring 
ferromagnctism but they can, in principle, take different values. We shall focus 
on the usual uniform case Jij — J, that is to say the clean model, and the 
case in which the J^'s are quenched random variables drawn from a probability 
distribution with positive support, a dirty case. The latter defines the random- 
bond Ising model (RBIM) . One can also add a random field term to the uniform 
model and thus construct the random field Ising model (RFIM). 

A totally random spin configuration in which Si = ±1 with probability a half 
is a realization of a site percolation configuration with p = 1/2 and, in d = 2, 
neither positive nor negative clusters percolate. The identification is achieved 
by defining site occupation variables rii = {si -I- l)/2 = 0, 1. 

The static properties of the Ising model (1) are more easily analyzed in the 
experimentally relevant canonical ensemble in which the system is coupled to a 
thermal reservoir at temperature T. A second-order continuous phase transition 
at a finite critical temperature, Tc, separates a high-temperature disordered 
paramagnetic phase from a low-temperature ferromagnetic ordered one. see 
Fig. 1. Tc depends on d and the statistics of the interaction strength. The 
nature of the two phases is understood since the development of the Curie- Weiss 
mean-field theory while the critical phenomenon has been accurately described 
with the help of the renormalization group, that allows for the computation of 
all critical exponents [4]. 

2.2.2 The Potts model 

A natural generalization of the Ising model is the Potts model [5] in which 
the spins cr^ take q integer values from 1 to q. The Hamiltonian is given by 



where the sum is over nearest-neighbours on the lattice. Either uniform, Jij = J 
(pure case), or bimodal P{Jij) = p5{Jij — Ji) -I- {l—p)5{Jij — J2) (random case) 
interactions are usually studied. In d = 2 the transition, discontinuous for g > 4 
and continuous for g < 4, occurs at (e'^^''^ — l){e^'''^'^ — 1) = g in the random 
case with p =1/2. Tc for the pure limit is recovered by setting Ji = J2. 

Soap films and general grain growth [6] are two physical systems mimicked 
by Potts models. In 2d these physical systems are made of polygonal- like cells 
separated by thin walls endowed with line energy; different 'colour' domains 
correspond to different cells. 

2.3 Microscopic dynamics 

Once the Hamiltonians have been defined, we need to assign updating rules 
to the microscopic variables. Classically, the spins do not have an intrinsic 
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Figure 1: An equilibrium free-energy landscape with a second-order phase tran- 
sition between two phases, a single-valued one with one minimum, and a bi- 
valued one with minima related by parity. Below arc three snapshots of a 
spontaneously broken state a,t g < critical configuration at g = gc, and a 
disordered one at g > g^, in the 2dlM. Black and white dots represent up and 
down spins, respectively. 

dynamics. Quite generally, we are interested in the evolution of physical systems 
in contact with an environment that provides thermal agitation and dissipation. 
Thus, the evolution is assumed to be stochastic [7] and it does not conserve the 
total energy. If we want to let the system reach equilibrium at sufficiently long 
times we use dynamic rules that satisfy detailed balance: 

WiC ^ C')PegiC) - W{C' ^ C)Peg(C") , (3) 

with W{C — !• C") the transition probability from configuration C to configura- 
tion C" and Peg(C) = e"'^^^'^'^^,, e-'^^C^") the Boltzmann weight. Finally, 
one has to decide whether there arc conserved quantities, for instance the global 
magnetization, M = X]i=i or a local one, e.g. mioc = Si + sj with i and j 
nearest-neighbours. The existence of conserved quantities puts constraints on 
the microscopic updates. In the context of Ising and Potts models two types of 
dynamics are common: 

Non-conserved order parameter, there are no constraints on the stochastic moves 
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and two examples are Monte Carlo or Glauber dynamics for Ising spins [8]. 
Conserved order parameter; it mimics particle-hole exchanges in lattice gas mod- 
els obtained from a mapping of the Ising model and the local magnetization is 
conserved. Exchanges of up and down spins are proposed and these are accepted 
with a probabilistic law (Kawasaki dynamics [8] ) . 

Kinetic constraints Spin updates are allowed (and realized with some stochastic 
rule) only when a chosen constrained is satisfied, typically by neighbor variables. 

Dynamical systems are defined by dynamic rules that do not necessarily 
satisfy detail balance and have no reference to any equilibrium energy (or free- 
energy) that they should minimize. Up to what extent the stochastic and deter- 
ministic evolution of some macroscopic systems are similar or even equivalent 
at some length and time scales is a relevant question. Some issues one would 
like to understand are under which conditions collective behaviour emerges in 
extended dynamical systems with short-range interactions and local chaotic dy- 
namics; and whether thermodynamic and statistical mechanics concepts apply 
in some dynamic regimes, in particular, whether an energy function and a tem- 
perature can be identified in deterministic dissipative systems. In spatially 
extended dynamical systems the effective noise strength, that may play the role 
of temperature, is internally generated rather than imposed by an environment. 
Some lattice models of coupled chaotic maps present dynamic phase transitions 
between dynamic phases that can be associated to a disorder state and a bi- 
valued dynamically ordered one by using a spin representation. This is the case 
of the Miller- Huse coupled map lattice [9], for example, the transition of which 
presents many similarities with critical phenomena. 

Problems that originate in other sciences, suchlike social studies, informatics, 
etc. can sometimes be formulated by using stochastic dynamic rules that do not 
satisfy detail balance. In some cases, collective behaviour gives rise to dynamic 
phase transitions between phases [3] and by choosing adequately the dynamic 
rules one may also encounter a dynamic transition similar to the ones mentioned 
above. This is the case, for instance, of some variants of voter models [10]. 

2.4 The spiral model 

Spin kinetically constrained models capture many features of real glass- 
forming systems (see [11] for reviews). These models display no thermodynamic 
singularity: their equilibrium measure is simply the Boltzmann factor of inde- 
pendent spins and correlations only reflect the hard core constraint. Bootstrap 
percolation arguments allowed C. Toninelli et al. to show that, when defined on 
finite dimensional lattices, these models do not even have a dynamic transition 
at a particle density that is less than unity in the thermodynamic limit [12]. 
On Bethe lattices instead a dynamical transition similar to the one predicted 
by the mode coupling theory might occur [13]. The quest to define a finite di- 
mensional kinetically constrained model with a (discontinuous) transition was 



6 



p=0.99 , 1=10* p=0.99 , l=10' 




Figure 2: Left: The neighbouring sites determining the frozen or free to move 
character of the center site Centre and right panel: Configuration of the 

system at two times {t = 10^ and t = 10^) in a quench to p = 0.99. Black and 
white sites are frozen particles and vacancies, respectively. Red and green sites 
are particles and holes that can be updated. Figure taken from [15]. 

positively answered by C. Toninelli et al who constructed the so-called spiral 
model, a finite dimensional kinetically constrained model with an ideal glass- 
jamming dynamic transition at a particle density that is different from one [14]. 
A binary variable = 0, 1 is defined on the sites of an L x L periodic 
bi-dimensional square lattice, riij = 1 if a particle occupies the site and = 
otherwise. One defines the couples of neighbouring sites (see Fig. 2-left): 
+ + l,j + 1), north east (NE) couple; 

(i + 1, j),(i + 1, j — 1), east south (ES) couple; 
— l),(t — 1, j — 1), south west (SW) couple; 

(i — — 1, j + 1), west north (WN) couple. 

Whether the variable may be updated or not depends on the configuration 
on these pairs: if the sites belonging to at least two consecutive couples (namely, 
NE+ES or ES+SW or WN+NW or WN+NE) are empty site (ij) can be up- 
dated (either emptied or filled). Otherwise it is blocked. Each site is coupled to a 
particle reservoir in such a way that particles can enter or leave the sample from 
its full volume. Defining My = if site is frozen and = 1 otherwise 
the updating rule can be expressed in terms of the transition rates W{nij\n'^j) 
to add or remove a particle as Wp(0|l) = My p and Wp(l|0) = Afy (1 —p). The 
model can be regarded as a two-level system described by = — J^ij (2ny — 1) 
at the inverse temperature (3 ~ (l/2)ln[p/(l — p)]. When p ^ 1/2 the in- 
verse temperature vanishes /? ^ whereas for p ^ 1 it diverges /3 ^ oo. The 
Bernoulli measure implies p = (riij) = p for the equilibrium density of particles. 
In equilibrium at zero temperature {p — > 1), the lattice is full while at infi- 
nite temperature (p — > 1/2) it is half-filled. The existence of a blocked cluster 
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at p > Pc < 1 was shown by studying the critical p above which an equihb- 
rium configuration cannot be emptied and this coincides with the threshold for 
2d directed percolation {pc ~ 0.705). At the transition an infinite cluster of 
blocked particles exists and the density of the frozen cluster is discontinuous at 
Pc- Many aspects of super-cooled liquids slowing down are reproduced close (but 
below) Pc [14]. In Sect. 7 we shall summarize results in [15] where we studied, 
in particular, the dynamics following a quench from pq < pc to p > p^. 

3 Geometric description of phase transitions 

Geometric domains are ensembles of connected nearest-neighbour sites with 
spins pointing in the same direction. Their area is the number of spins belonging 
to the domain. By lumping together with a certain temperature-dependent 
probability neighboring spins in the same spin state, spin models can be mapped 
onto percolation theory. The resulting Fortuin-Kasteleyn spin clusters [16] built 
using — 1 — e^^'^'J percolate at the critical temperature, and their percolation 
exponents coincide with the thermal ones. In this way, a purely geometrical 
description of the phase transition is achieved. In the rest of this Section we 
sum up the definitions of the main geometric objects used in such a static 
description of critical phenomena. 

3.1 Definitions 

A number of linear dimensions of a cluster can be defined. The radius of 
gyration is 

Rl-S-'Y.\'-~rcra\' (4) 
i=l 

with the centre of mass position given by rem — s~'^ X^i^i 

Take the cluster site with the largest (smallest) y component and among 
these the one with the largest (smallest) x component. These sites are the 
two 'end-points' of the cluster. The end-to-end length of the cluster, i, is the 
Euclidean distance between these two points. 

There are also several ways of defining the surface of a cluster [2] . Different 
definitions are relevant to different applications that have to do, for example, 
with the adsorption of particles on the surface. Without entering into all the 
zoology let us recall some of these definitions. 

The external border of a cluster of occupied sites is the set of vacant sites 
that are nearest-neighbours to sites on the cluster and arc connected to the 
exterior by (a) nearest and ncxt-to-nearest neighbours or (b) just next nearest 
neighbours. Figure 3 (a) and (b) show the external borders of a cluster of 
occupied sites thus defined. Note that in (a) there arc two 'inner' sites that 
belong to the external border and that are connected to the outside by a narrow 
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neck of width \f2a with a the lattice spacing while in (b) these two sites do not 
belong to the border and the 'fjord' has been excluded. 

The hull is the envelop of a cluster, meaning the ensemble of cluster sites 
obtained by using a biased walk that encircles the cluster on its left and on its 
right linking the lowest lying site to the left and the highest lying site to the 
right, see Fig. 4. The hull of a cluster of occupied (vacant) sites lies on occupied 
(vacant) sites. 

The hull-enclosed area is the total area within the hull (including the sites 
on the hull). It therefore ignores the type of site (occupied or vacant, spin up 
or spin down) that lies within and counts them all on equal footing. 

3.2 Number densities 

We call n(s) the number of objects (hulls, clusters, etc.) with s sites per 
unit number of lattice sites. Scaling theory, some exactly solvable cases (^e.g. 
the Bcthe lattice), renormalization group arguments [2], Coulomb gas mappings, 
conformal field theory calculations, stochastic Loewner evolution techniques [17, 
18] and numerical experiments [19] suggest 

ri(s) ^ s^'^ f[Os'^] for large s , (5) 

with T and a two control parameter and lattice-independent exponents that do 
depend on the space dimensionality and 9 measuring the distance from critical- 
ity, say 6 = \g — gd with g the control parameter. The function f(z) approaches 
a constant for |z| ^ 1 and it falls-off rapidly for |2| ^ 1; it thus provides a 
cut-off with a unique cross-over size = [5 — ffd"^^'^ such that clusters with 
s < are effectively critical while those with s > are not. At criticality 
{g = gc) the factor f[9s'^] = /[O] that suppresses large objects is absent and 
clusters of all sizes arc present. 

The exponents a and r determine all critical exponents of percolation through 
scaling relations. As regards thermal phase transitions, once the relevant Fortuin- 
Kasteleyn clusters arc constructed and analysed their two independent expo- 
nents a and r determine the entire set of thermal critical exponents (see, how- 
ever, [20]). 

3.3 Fractal exponents 

Geometric objects at criticality have a fractal behaviour, that is to say [2] 

s oc £^ . (6) 

s is the mass of the object (cluster area, its external border, its hull), ^ is 
the linear size of the cluster, say its radius of gyration, and D is the fractal 
dimension. The fractal dimension D is related to the exponent r in eq. (5) via 
T — d/ D + 1 where d is the space dimension. 
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Figure 3: Square bi-dimcnsional lattice. Filled (empty) circles represent occu- 
pied (empty) sites. Red crosses indicate external border of the cluster of filled 
sites including (a), or nor (b), vacant sites that are connected to the exterior by 
a next nearest neighbour (a 'diagonal'). 



The fractal exponents depend strongly on the precise definition of the ge- 
ometric object. For instance, Salcur and Duplanticr showed that 2d critical 
percolation clusters have De = 4/3 and Dh = 7/4 where H stands for hull and 
E for external border connected to the exterior by nearest-neighbour vacant 
sites only [21] (case (b) in Fig. 3). 

In [22] we studied the relation between areas and perimeters 

A cx (7) 

with p defined as the number of broken bonds on the external border or, equiva- 
lently, with a small variant of the Grossman- Aharony construction [2] including 
fjords. Note that ii A (x and p oc i^^ with £ the linear size of the object, 
say the radius of gyration, then a = Da/ Dp. 

4 Coarsening 

Take a system in equilibrium in the symmetric (positive g ^ gc, high tem- 
perature) phase and quench it into the symmetry breaking (negative g — gc, 
low temperature) phase through a second order phase transition. For concrete- 
ness, we focus on problems with two equilibrium states related by Z2 symmetry. 
Once set into the ordered phase the system locally selects one (among the two 
possible) equilibrium configurations. However, different 'states' are picked up at 
different locations and topological defects in the form of domain walls are cre- 
ated. In the course of time the patches of ordered regions tend to grow while the 
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Figure 4: Same eluster as in Fig. 3. Eneircled are the sites on the cluster that 
belong to its hull. The hull-enclosed area is everything that lies within the hull. 

density of topological defects diminishes. In an infinite system this coarsening 
process goes on forever. 

An example of the above process is given by an Ising magnet. An initial 
condition in equilibrium at, for instance, infinite temperature is just a random 
configuration in which each spin takes one of its possible values with probability 
1/2. After a quench below the critical temperature, T < Tc, the ferromagnetic 
interactions tend to align the neighbouring spins in 'parallel' direction (same 
value of the spin) and in the course of time domains of the ordered phases form 
and grow. At any finite time the configuration is such that both types of domain 
exist. Under more careful examination one reckons that there are some spins 
reversed within the domains. These 'errors' are due to thermal fluctuations 
and are responsible of the fact that the magnetization of a given configuration 
within the domains is smaller than one and close to the equilibrium value at 
the working temperature (apart from fluctuations due to the flnite size of the 
domains) . At each instant there are as many spins of each type (up to fluctuating 
time-dependent corrections that vanish in the infinite size limit). As time passes 
the typical size of the domains increases in a way that we shall discuss below. 

Coarsening from an initial condition that is not correlated with the equi- 
librium state and with no bias field does not take the system to equilibrium 
in finite times with respect to a function of the system's linear size, L. More 
explicitly, if the growth law is a power law [see eq. (10)] one needs times of the 
order of L^'' to grow a domain of the size of the system. This gives a rough 
estimate of the time needed to take the system to one of the two equilibrium 
states. For any shorter time, domains of the two types exist and the system is 
out of equilibrium. 

We thus wish to distinguish the relaxation time, tr, defined as the time 
needed for a given initial condition to reach equilibrium with the environment. 
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Figure 5: Snapshot of the 2d Ising model at a number of Monte Carlo steps 
after a quench from infinite to a sub-critical temperature. Left: the up and 
down spins on the square lattice are represented with black and white sites. 
Right: the domain walls are shown in black. Figure taken from [22]. 

from the decorrelation time, td, defined as the time needed for a given config- 
uration to decorrelate from itself. At T < Tc the relaxation time of any initial 
condition that is not correlated with the equilibrium state diverges with the 
linear size of the system. The self-correlation of such an initial state evolving at 
T < Tc decays as a power law and although one cannot associate to it a decay 
time as one docs to an exponential, one can still define a characteristic time 
that turns out to be related to the age of the system. 

In contrast, the relaxation time of an equilibrium magnetized configuration 
at temperature T vanishes since the system is already in equilibrium while the 
decorrelation time is finite and given by ~ |T — Td^"^"''. 

The lesson to learn from this comparison is that the relaxation time and the 
decorrelation time not only depend on the working temperature but they also 
depend strongly on the initial condition. Moreover, the relaxation time depends 
on (N, T) while the decorrelation time depends on (T, tw)- For a random initial 
condition one has 

r finite T>Tc, 
[ L^" T <Tc. 

4.1 Critical coarsening 

Right after a quench to the critical point the system starts to evolve out of 
equilibrium towards a target equilibrium configuration of the kind shown in the 
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central snapshot in Fig. 1. The non-equihbrium relaxation shows coarsening 
features with the growth of ordered structures with a typical linear length given 
by Rc(t) ~ t^/'''"' , and features of equilibrium dynamics, as the fact that Zeq is 
the equilibrium dynamics exponent. Interestingly enough, the early relaxation 
exhibits universal scaling properties characterized by usual static exponents as 
well as the dynamic one. This and other details of the critical equilibrium 
and out of equilibrium dynamics can be computed with renormalization group 
techniques [23, 24]. 

4.2 Dynamic scaling hypothesis 

The dynamic scaling hypothesis states that at late times and in the scaling 
limit [1] 

r » , R{9, t) » ^{g) , r/R{g, t) arbitrary , (8) 

where r is the distance between two points in the sample, r = \x — x'\, and 
S,{g) is the equilibrium correlation length that depends on all parameters (T and 
possibly others) collected in g, there exists a single characteristic length, R{g, t), 
such that the domain structure is, in statistical sense, independent of time when 
lengths are scaled by R{g, t). Time, denoted by t, is typically measured from the 
instant when the critical point is crossed. In the following we ease the notation 
and write only the time-dependence in R. This hypothesis has been proved 
analytically in very simple models only, such as the one dimensional Ising chain 
with Glauber dynamics or the Langevin dynamics of the d-dimensional 0{N) 
model in the large N limit. In the vast majority of coarsening systems the 
dynamic scaling hypothesis applies. Still, a few counter-examples where two 
lengths grow in competition arc also known [1] . 

In practice, the dynamic scaling hypothesis implies that the real-space cor- 
relation function should behave as 

C(r,t) EE 1^,, ~ F[r/Rit)] . (9) 

See Fig. 6 for a test of the scaling hypothesis in the MC dynamics of the 2d Ising 
and Potts model (see the caption for details). In order to fully characterise the 
correlation functions one then has to determine the typical growing length, R, 
and the scaling functions, fc, F, etc. 

4.2.1 The grov^ring length 

It turns out that R can be determined with semi-analytic arguments and 
the predictions arc well verified numerically - at least for clean system. In pure 
and isotropic systems the growth law is 

Rit) = A i^/^" (10) 
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Figure 6: Check of dynamic (super-)scaling in the space-time correlation of 
the Potts model quenched from Tq ^ oo to Tf = Tc{q)/2. Data taken at 
several times (from t = 2^ MCs to 2^^ MCs) for g = 2. 3 and 8, with and 
without weak disorder are shown. Distance is rescaled by the R{t) obtained 
from C{R,t) ~ 0.3. Inset: R^{t) against t in a double logarithmic scale for 
g = 2, 3 and 8 (from top to bottom). The characteristic length, related to 
the average domain radius, depends weakly on q and T for the pure model 
(through the pre- factor) . Within the time window explored there are still some 
small deviations from the R{t) ~ t^/^ expected law in the case q = 8. When 
weak disorder is introduced (not shown), the growth rate is greatly reduced and 
deviates from the linear behavior, see Sect. 5.2. Figure taken from [25]. 

with Zd the dynamic exponent (see [1]) and A a material/model dependent pref- 
actor that weakly depends on temperature and other parameters. In curvature 
driven Ising or Potts cases with non-conserved order parameter the domain walls 
have finite width and = 2. For systems with continuous variables such as 
rotors or XY models and non-conserved order-parameter dynamics, a number 
of computer simulations have shown that domain walls are thicker and Zd = 4. 
The effects of temperature enter only in the parameter A and, for clean systems 
of Ising type, growth is slowed down by temperature since thermal fluctuation 
tend to roughen the interfaces thus opposing the curvature driven mechanism. 
Let us list some special cases below and sketch how Zd can be estimated. 
Clean one dimensional cases with non-conserved order parameter 

In one dimension, a space-time graph allows one to view coarsening as the 
diffusion and annihilation upon collision of point-like particles that represent 
the domain walls. In the Glauber Ising chain with non-conserved dynamics one 
finds that the typical domain length grows as t^^^ while in the continuous case 



14 



the growth is only logarithmic, Int. 
Non-conserved curvature driven dynamics (d> 2) 

The time-dependent Ginzburg-Landau model allows us to gain some insight 
on the mechanism driving the domain growth and the direct computation of 
the averaged domain length. This is a stochastic partial differential equation 
on a coarse-grained order parameter field with a deterministic force that is phe- 
nomenologically proposed to derive from a Ginzburg-Landau type free-energy. 
In clean systems temperature does not play a very important role in the domain- 
growth process, it just adds some thermal fluctuations within the domains, as 
long as it is smaller than Tc- In dirty cases instead temperature triggers thermal 
activation. 

We focus first on zero temperature clean cases with Ising-like symmetry. At 
T = the GL equation is just a gradient descent in a (free-)energy landscape, 
F. Two terms contribute to F: an elastic energy (Vt/))^ which is minimized 
by flat walls if present and a bulk-energy term that is minimized by constant 
field configurations, say — ±00 in an Ising-likc case. If the walls are sharp 
enough, that is to say their width remains finite when the distance between them 
diverges, interface-interface interactions can be neglected and the minimization 
process implies that regions of constant field grow and get separated by flatter 
and flatter walls. Within this scenario one can easily derive the Allen-Cahn 
equation [26] that states that the local wall velocity is proportional to the local 
geodesic curvature and is normal to the wall pointing in the direction of reducing 
curvature: ^ 

H = dtn\d, = —V ■ fi fi = — — nn , (11) 

ZTT 

in all d. This equation yields an intuition on the typical growth law in such 
processes. Take a spherical wall in any dimension. The local curvature is k = 
{d— l)/R where R is the radius of the sphere within the wall. Equation (11) is 
recast as dR/dt = -X{d - 1)/R that implies R'^{t) = i?2(0) - 2X{d - l)t and R 
decreases as t^^^. This calculation implies that all hull-enclosed areas decrease 
in time with the same law, independently of their own size and of all others. 
Interestingly enough, similar results are obtained for the area loss of a square 
domain embedded in a sea of the opposite sign using Monte Carlo or Glauber 
dynamics at zero temperature for a system defined on a square lattice [27]. 

The above results do not imply that all domains shrink as well since some 
gain size from the disappearance of inner objects. Temperature effects are sim- 
ply taken into account by a material and T-dependent proportionality constant 
A in the Allen-Cahn equation. 

In d = 2 the time dependence of the area contained within any finite hull 
on a flat surface is derived by integrating the velocity around the envelope and 
using the Gauss-Bonnet theorem: 
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where 9i are the turning angles of the tangent vector to the surface at the n 
possible vertexes or triple junctions between domains of different colour (we are 
now generalizing the discussion to models with multi-valued equilibrium states 
suchlike the Potts model with q > 2). In the Ising q ^ 2 model, J^i^i — 
since there are no such vertexes and we obtain dA/dt = — A for all hull-enclosed 
areas, irrespective of their size. If, instead, like in soap froths, there is a finite 
number of such vertexes with an angle of 27r/3, that is, 9i = ir/S^Wi (for highly 
anisotropic systems, as the Potts model, the angles are different from 2tt/3). The 
above equation thus reduces to the von Neumann law [28] for the hull-enclosed 
area An with a hull with n-turning angles: 

(in the Ising case n = 0). Whether a cell grows, shrinks or remains with constant 
area depends on its number of sides being, respectively, larger than, smaller than 
or equal to 6. 

Conserved order parameter: bulk diffusion 

A different type of dynamics occurs in the case of phase separation (a water 
and oil mixture ignoring hydrodynamic interactions or a binary alloy). In this 
case, the material is locally conserved, i.e. water does not transform into oil 
but they just separate. The main mechanism for the evolution is diffusion of 
material through the bulk of the opposite phase. After some discussion, it was 
established, as late as in the early 90s, that for scalar systems with conserved 
order parameter 2d = 3 (see [29, 8]). 
Role of disorder: thermal activation 

The situation becomes less clear when there is quenched disorder in the form 
of non-magnetic impurities in a magnetic sample, lattice dislocations, residual 
stress, etc. Qualitatively, the dynamics are expected to be slower than in the 
pure cases since disorder pins the interfaces. In general, based on an argument 
due to Larkin (and in different form to Imry-Ma) one expects that in d < 4 
the late epochs and large scale evolution be no longer curvature driven but 
controlled by disorder. 

The argument to estimate the growth law in dirty systems is the following. 
Take a system in one equilibrium state with a domain of linear size R of the 
opposite equilibrium state within it. This configuration could be the one of an 
excited state with respect to the fully ordered one with absolute minimum free- 
energy. Call AF{R) the free-energy barrier between the excited and equilibrium 
states. The thermal activation argument yields the activation time scale for the 
decay of the excited state (i.e. erasing the domain wall) 

tA - TO eAF(i?)/(/c,T) _ (14) 

For a barrier growing as a power of i?, AF{R) - T(T, J)R'^' (where J represents 
the disorder) one inverts (14) to find the linear size of the domains still existing 
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at time t, that is to say, the growth law [30] 



ksT , t 



) 



i/V 



T(r,j) To 



(15) 



All smaller fluctuation would have disappeared at t while typically one would 
find objects of this size. The exponent ip is expected to depend on the dimen- 
sionality of space but not on temperature. In 'normal' systems "0 should be 
d — 1 - the surface of the domain. The pre-factor T is expected to be weakly 
temperature dependent. 

To extend this result to the actual out of equilibrium coarsening situa- 
tion one assumes that the same argument applies out of equilibrium to the 
re-conformations of a portion of any domain wall or interface where R is the 
observation scale. 

However, not even for the (relatively easy) random fcrromagnet there is 
consensus on the actual growth law [31]. We shall discuss a possible way out in 
Sect 5.1. In the case of spin-glasses, if the mean-field picture with a large number 
of equilibrium states were realized in finite dimensional models, the dynamics 
would be one in which all these states grow in competition. If, instead, the 
phcnomenological droplet model applied, there would be two types of domains 
growing, R{t) ^ {hity^'^ and the dimension of the bulk of these domains should 
be compact but their surface should be rough with fractal dimension ds > 
d-1 [32]. 

4.2.2 Scaling functions 

The scaling functions are harder to obtain. Indeed, there is no system- 
atic method to derive them and all approximations dealt with are not accurate 
enough. For a much more detailed discussion of these methods see the review 
articles in [1]. Still, the super- universality property seems to be correct as we 
discuss in Sect. 5.2. 

5 Dynamics of (weakly) random systems 

In this Section I explain that the existence of a static typical length yields 
a natural crossover between the clean growth law, say R ~ [A(T)t]^/^, and the 
activation-ruled one, R ~ [fcsT/T Ini]!/'^, and how these two regimes might not 
be sufficiently separated in numerical and experimental measurements, being 
easily confused with a disorder and temperature dependent power [33] . 

I also discuss the super-universality hypothesis [32] and list some of its 
checks. 

5.1 Crossover in the growing length 

Numerical simulations of dirty systems tend to indicate that the growing 
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length is a power law with a disorder and T-dependent exponent. This can be 
due to the effect of a disorder and T-dependent cross-over length, as explained 
in [33]. For concreteness, let us assume that the width of the quenched random 
distribution is characterised by a parameter J and that the cross-over length 
depends on J/T. We call the latter Lt by absorbing the J dependence in T. 
The proposal is that below Lt the growth process is as in the clean limit while 
above Lt quenched disorder is felt and the dynamics are thermally activated 
above barriers that are usually taken to grow as a power of the size leading to 
eq. (15): 

/ [KtW^- for i?W«LT, 

^ ^ ^ I [fci5r/T(r)lnt]i/^ for i?(t)»LT. ^ ^ 

These growth-laws can be first inverted to get the time needed to grow a given 
length and then combined into a single expression that interpolates between the 
two regimes: 

t{R) - e'^R/LT)^R-^ (17) 

where the relevant T-dependent length-scale Lt has been introduced. Now, by 
simply setting t{R) ~ i?^*^-^^ one finds 

z(T) - ~ z/L% for times such that f^''-" \nt ~ ct . (18) 

Similarly, by equating t{R) - exp(i?^'(^) /keT) one finds that V'(r) is a decreas- 
ing function of T approaching ip at high T . 



5.2 Super-universality 

In its first version the super-universality hypothesis states that in cases in 
which temperature and quenched disorder arc 'irrelevant' in the sense that they 
do not modify the nature of the low-temperature phase {e.g. it remains fer- 
romagnetic in the case of ferromagnetic Ising models with quenched random 
interactions) the scaling functions should not be modified [32]. Only the grow- 
ing length changes from the, say, curvature driven t^/^ law to an asymptotically 
slower law due to domain wall pinning by impurities. Tests of the disorder 
independence of the scaling functions of several correlations in quenches from 
T > Tc into the ordered phase in the M RFIM [34] and the 2d RBIM [35, 36] 
(perhaps excluding the possibility of having zero bonds, see the discussion in 
Henkel & Pleimling) give support to this hypothesis. 

We recently investigated two aspects of the super-universality hypothesis. 
On the one hand, we focused on the dependence on the initial condition of the 
space-time correlation scaling functions in the low temperature phase. With this 
aim, we studied the dynamics of the Potts model with different values of q after 
quenches from To ^ oo and Tq = Tc [25]. We found that the scaling functions 
are fully determined by the type of correlations present in the initial conditions 
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and basically of two types. All 2 < q < 4 clean and weak disordered Potts 
models quenched from equilibrium at To — > oo and clean g > 4 Potts models 
quenched from equilibrium at Tc (the transition is of first order) share their 
scaling function (see Fig. 6). Different functions are obtained after quenches 
from Tc in clean systems with 2 < q < 4 and dirty systems with q > 4; in 
all these cases the transition is of second order and the critical exponents, in 
particular 77, depend on q (although weakly) [37]. 

On the other hand, we analysed the geometric properties of areas and perime- 
ters in the 2c? RBIM and we discuss these results in the next section. 

6 Statistics and geometry of coarsening 

In [38, 22] we studied the distribution of domain areas, areas enclosed by do- 
main boundaries, and perimeters for curvature-driven two-dimensional Ising-like 
coarsening, employing a combination of exact analysis and numerical studies, 
for various initial conditions. We showed that the number of hulls per unit area, 
nh{A, t) dA, with enclosed area in the interval {A, A + dA), is described, for a 
disordered initial condition, by the scaling function 

nh{A,t)=2ch/{A + \htf , (19) 

where cn = l/87r^/3 « 0.023 is a universal constant and \h is a material pa- 
rameter. For a critical initial condition, the same form is obtained, with the 
same A^i but with Ch replaced by Ch/2. For the distribution of domain areas, we 
argued that the corresponding scaling function have the form 

nd(A, t) = {2)cd{\dtY-^/{A + Xdty , (20) 

where Cd and Ad are numerically very close to Ch and Xt respectively, and the 
exponent r is the one characterising the distribution of initial structures, crit- 
ical percolation or critical Ising. These results were extended to describe the 
number density of the length of hulls and domain walls surrounding connected 
clusters of aligned spins. These predictions were supported by extensive numer- 
ical simulations. We also studied numerically the geometric properties of the 
boundaries and areas. 

The derivation of eq. (19) is very easy. The number density of hull enclosed 
or domain areas at time t as a, function of their initial distribution is 

/•OO 

n(A,t)= / dA.,5{A~ A{t,A,))n{A,,U) (21) 
Jo 

with Ai the initial area and n{Ai^ti) their number distribution at the initial 
time ti. A{t,Ai) is the hull enclosed/domain area, at time t, having started 
from an area Ai at time ti. The number density of hull-enclosed areas in critical 
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Figure 7: Left: number density of hull-enclosed areas at different times after 
a quench from Tg ^ oo to T = in the 2dlM; numerical results are shown 
with points and the analytic prediction with lines. There is only one fitting 
parameter, A, that has been independently determined by a study of the space- 
time correlation. Right: comparison between the scaled number-density after a 
quench from Tq — > oo (above) and Tq = Tc (below). Figure taken from [22]. 

percolation and critical Ising conditions was obtained by Cardy and Ziff [39] 



These results are valid for Aq A L^, with Aq a microscopic area and the 
system size. Note also that we arc taking an extra factor 2 arising from the fact 
that there arc two types of hull enclosed areas, corresponding to the two phases, 
while the Cardy-Ziff result accounts only for clusters of occupied sites (and not 
clusters of unoccupied sites). nh{A,0) dA is the number density of hulls per 
unit area with enclosed area in the interval (A, A + dA) (we keep the notation 
to be used later and set t = 0). The adimensional constant Ch is a universal 
quantity that takes a very small value: Ch = l/{8^T^/3) « 0.022972. The time 
dependent area is given by the AUen-Cahn result, A{t) = Ai — \h{t — ti). By 
simple integration of eq. (21) one obtains eq. (19). 

It is interesting to note that the infinite temperature initial condition is not 
critical percolation on the square lattice. Still, it is very close to it an after 
coarse-graining p =1/2 becomes critical continuous percolation. The dynamics 
of the discrete model, thus, very quickly settle into critical percolation 'initial' 
conditions and this determines the distribution of large structures dynamically. 
This observation was used by Barros et al to interpret freezing of 2d Ising models 
at very low temperatures [40] . 

The calculation of the number density of domain areas cannot be done ana- 
lytically but a mean-field-like approximation that uses an expansion in powers 




critical percolation, 
critical Ising. 



(22) 
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of Ch yields eq. (20), a result that is very convincing since it agrees well with 
numerical simulations [38, 22]. 

The discussion on the geometric properties of critical objects suggests to 
study the fractal properties of growing structures. In [38, 22] we pursued this 
line of research by analysing the relation between areas and perimeters during 
coarsening. We summarize the results below. 

Several extensions of these results are: to the Ising model with conserved 
order parameter [41], the random bond Ising model [36], the Potts model [25], 
and coupled map lattices [42]. 



6.1 General picture 

The summary of results in, and picture emerging from, the study of coars- 
ening in 2d systems with Ising symmetry is the following. 

(i) We proved scaling for the number density of hull-enclosed areas in bi- 
dimensional zero-temperature curvature-driven coarsening in the Ising univer- 
sality class [38]. 

(ii) We argued that in these systems temperature effects are two-fold: they 
renormalize the pre-factor in the growth law and they introduce thermal do- 
mains that are distributed as in equilibrium. Their contribution to the total 
distribution of structures can be safely subtracted [22] . 

(iii) We obtained approximate expressions for the number density of domain 
areas and perimeters in curvature-driven coarsening in the 2d Ising universality 
class by using a mean- field-like analytic argument [22]. 

(iv) Wc checked results (i)-(iii) experimentally by analyzing the dynamics of 
a 2d liquid crystal [43]. 

(v) We derived the scaling functions of the number density of hull-enclosed 
areas, domain areas and perimeters for systems in the universality class of the 2d 
Ising model class with conserved order-parameter dynamics quenched from high- 
temperature by assuming that small structures behave roughly independently 
of each other [41]. 

(vi) Wc analysed the area number densities in the RBIM and we found 
that super-scaling holds once the growing length is modified to incorporate the 
slowing down due to pinning [36]. 

All these results are described by the following conjecture for the number 
density of domains and hull-enclosed areas [41] 

1/2 



A 

Rm 



1 



^ tW2^ (2r+l)/., 



(23) 



K^t) 



(this expression is exact for hull-enclosed areas in curvature-driven coarsening) . 
Zd is the dynamic exponent (in the RBIM case, it is the effective temperature- 
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dependent one, Zrf). t characterises the initial distribution. A way to derive 
eq. (23) is to assume that each time-dependent area is independently linked to 
its initial value by 

A^'-'/^it) ^ Ai'"^^ R{t) , (24) 

In all these cases the expressions that we obtained have two distinct limiting 
regimes. 

(vii) For areas that arc much smaller than the characteristic area, (t) , the 
distributions 'feel' the microscopic dynamics and thermal agitation while for 
areas that arc much larger than (t) the distributions are simply those of the 
initial condition. For conserved order parameter dynamics the Lifshitz-Slyozov- 
Wagner [44] behaviour is recovered after a quench from Tq — > oo and evolving at 
sufficiently low T. For critical Ising initial conditions the distribution of small 
areas does not approach the Lifshitz-Slyozov- Wagner. We conjectured that the 
reason is that our starting assumption, independence of domain wall motion for 
small domains, is not valid due to strong correlations in this case [41]. 

(viii) The distribution of the time-dependent areas that are larger than R^ (t) 
are, in all cases, the ones of critical continuous percolation (for all initial condi- 
tions in equilibrium at Tq > Tc) and critical Ising (Tq = Tc). 

We also studied the geometric properties and fractal dimensions of different 
objects during coarsening. We found that 

(ix) Small structures are compact and tend to have smooth boundaries, that 
are pretty close to circular (A ~ p^) in the conserved order-parameter case and 
a little bit less so for non-conserved order-parameter dynamics. 

(x) The long interfaces retain the fractal geometry imposed by the equilib- 
rium initial condition. 

In [43] we performed experiments in a 2d liquid crystal that evolves through 
the formation of domains of two chiralities. The analysis of data confirmed 
that the dynamics is of curvature-driven type, in the sense that the correlation 
functions scale with R ^ t^l"^ and, moreover, after a careful analysis of noise 
originated in the data acquisition procedure, we obtained a hull-enclosed area 
distribution function that agrees well with the theoretical prediction. 

The geometrical analysis of coarsening cells in the 2d Potts model will be 
presented in [25]. One might expect that other dynamic systems, once in the 
Ising non-conserved order parameter universality class, such as some voter mod- 
els [10], should have similar area distributions. The case of a 2d Edwards- 
Anderson model has not been analyzed yet. 

As future work we plan to investigate how these results extend to 3d systems 
by analyzing volume distributions as well distributions on 2d cuts. Three 

dimensional phase separation of a binary mixture is realized by different physical 
systems, in particular the one studied in [45], that promises to be a nice material 
where to confront the results of analytics and simulations to experiments. 
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7 Coarsening in the spiral model 



Among other features, in [15] we studied the dynamics of the spiral model 
after a sudden quench, in which a completely empty configuration is evolved 
from time t = onwards with the dynamic rule specified by a different value 
of p > pq. This procedure is similar to the temperature quench of a liquid. 
For p < pc, after a non equilibrium transient the system attains a non-blocked 
equilibrium state. The relaxation time for attaining such a state diverges as 
p ^ Pc- In the critical case with p = Pc, the system approaches the blocked 
equilibrium state by means of an aging dynamics similar to that observed in 
critical quenches of ferromagnetic models. Freezing is not observed because the 
dynamic density p{t) is always smaller than the critical one pc = Pc, at any 
finite time. Something different happens for filling at p > Pc- The system keeps 
evolving but a blocked state is never observed despite the fact that the density 
of particles exceeds Pc = Pc at long enough times. This is no surprise since the 
dynamics are fully reversible: for any p, each configuration reached dynamically 
can always evolve back to the initial empty state by the time reversed process, 
although with a very low probability. Therefore, a blocked state cannot be 
dynamically connected to the initial empty state. The dynamics of the spiral 
model at p > Pc strongly resembles coarsening in ferromagnets, as illustrated 
in the centre and right panels in Fig. 2. In the p ~ 1 limit in which the 
dynamics can be analyzed scmi-analytically the system coarsens by forming 
longer and longer onc-dimensional objects of vacancies that are almost frozen 
but not completely since they could be destroyed by the boundaries. The density 
of these objects is irrelevant in the large size limit and the density of particles 
can asymptotically approach 1 although the system is never blocked. 
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